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ABSTRACT 

There is currently much interest in the possible presence of intermediate-mass black 
holes in the cores of globular clusters. Based on theoretical arguments and simulation 
results it has previously been suggested that a large core radius - or particularly a large 
ratio of the core radius to half-mass radius - is a promising indicator for finding such 
a black hole in a star cluster. In this study ./V-body models of 100 000 stars with and 
without primordial binaries are used to investigate the long-term structural evolution 
of star clusters. Importantly, the simulation data is analysed using the same processes 
by which structural parameters are extracted from observed star clusters. This gives a 
ratio of the core and half-mass (or half-light) radii that is directly comparable to the 
Galactic globular cluster sample. As a result, it is shown that the ratios observed for 
the bulk of this sample can be explained without the need for an intermediate-mass 
black hole. Furthermore, it is possible that clusters with large core to half-light radius 
ratios harbour a black-hole binary (comprised of stellar mass black holes) rather than 
a single massive black hole. This work does not rule out the existence of intermediate- 
mass black holes in the cores of at least some star clusters. 

Key words: stellar dynamics — methods: N-body simulations — stars: evolution — 
binaries: close — globular clusters: general — open clusters and associations: general 



1 INTRODUCTION 

The situation regarding the growing body of evidence 
1 that some globular clusters (GCs) may be harbouring 
intermediate-mass black holes (IMBHs) has been summa- 
rized recently by Baumgardt, Makino & Hut (2005). This 
evidence includes taking the relationship found between the 
masses of supermassive black holes (BHs) and the bulge 
masses of the host galaxies (Magorrian et al. 1998) and ex- 
trapolating to globular cluster masses (Kormendy & Rich- 
stone 1995). For a typical globular cluster, such as M15 (van 
der Marel 2001), this gives a BH mass of ~ 10 3 Mq. Sitting 
conveniently between the supermassive and stellar-mass BH 
regimes - where the latter includes BHs of ~ 50 Mq or less - 
the IMBH tag arises naturally. The existence of such BHs is 
backed up by the iV-body simulations of Portegies Zwart et 
al. (2004) showing that possible progenitors (main-sequence 
stars of ~ 10 3 Mq) can be created through runaway merg- 
ers of massive stars in young clusters. Detection is possible 
through the measurement of central velocity dispersions in 
globular clusters but this is a challenging process (Baum- 
gardt, Makino & Hut 2005; Trenti 2006). To date this has 
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led to suggestions of an IMBH in the core of M15 (Gerssen 
et al. 2002) and in the core of Gl (Gebhardt, Rich & Ho 
2002). However, Baumgardt et al. (2003a, 2003b) subse- 
quently used iV-body simulations to show that the inferred 
non-luminous central mass could instead be a central con- 
centration of stellar-mass BHs, white dwarfs and neutron 
stars (but see also Gebhardt, Rich & Ho 2005). 

Notwithstanding the lack of direct confirmation that 
IMBHs do reside in the cores of GCs, study into the ramifica- 
tions of such a scenario has progressed. Importantly, Baum- 
gardt, Makino & Hut (2005) have shown that a GC with an 
IMBH in the core will be observed to have a relatively fiat 
central surface brightness profile and consequently a larger 
measured core-radius compared to a GC without an IMBH. 
This result has been followed up by Trenti (2006) who sug- 
gests that the ratio of the core radius, r c , to the half-mass 
radius, r^, of a dynamically-evolved cluster can be used to 
infer the presence of an IMBH. Trenti (2006) combines re- 
sults from a variety of TV-body simulations (Heggie, Trenti 
& Hut 2006; Trenti, Heggie & Hut 2007; Trenti et al. 2007). 
These show that r c /rh ~ 0.02 for clusters composed initially 
of single stars only, r c /r^ ~ 0.05 for clusters with primor- 
dial binaries, and r c /r^ ~ 0.3 for clusters with an IMBH. 
These values are taken when the model clusters are relaxed 
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systems and the core-collapse phase has ended. In compar- 
ison, observations of Galactic GCs show a distribution of 
r c /rh extending from 0.1 — 1.0 with a peak at about 0.5 
(Fregeau et al. 2003). From a theoretical viewpoint Heggie 
et al. (2006) examine how the r c /rh ratio varies with the 
BH mass. This also suggests that a star cluster observed to 
have a large core radius presents the most promising target 
for finding an IMBH, in the sense that large mass implies 
large core radius. As with the above results this argument 
is only valid in the post-collapse regime. 

A recurring issue with A-body simulations of star clus- 
ter evolution is that the models are generally idealized in 
some way (or ways) that prohibits direct comparison to real 
clusters. The simulations of Heggie, Trenti & Hut (2006), 
Trenti, Heggie & Hut (2007) and Trenti et al. (2007) were 
restricted to initial particle numbers of No = 20 000 or less 
and assumed equal-mass stars. As pointed out by Trenti 
(2006) these results can be scaled to GC particle numbers 
(No ~ 10 5 — 10 6 ) but only by also neglecting stellar evolu- 
tion. Simulations performed by Baumgardt & Makino (2003) 
and Baumgardt, Makino & Ebisuzaki (2004) included parti- 
cle numbers up to 131 072 stars, a mass spectrum and stellar 
evolution. However, primordial binaries were not included. 
Another key factor is that one must be sure to compare like- 
with-like when using model and real data. Specifically this 
relates to use of the core radius, half-mass (or half-light) 
radius, and the half-mass relaxation timescale, t T h- 

Considering the growing interest in IMBHs it is only 
natural that attempts are being made to isolate key ob- 
servational tests for their existence. Unfortunately, in this 
paper, it is shown that r c /r^ cannot readily be used as such 
a test. This is based on a series of A-body models of 100 000 
stars with and without primordial binaries. The models in- 
clude a full mass spectrum, stellar and binary evolution, 
and account for the tidal field of the Galaxy. The models do 
not include IMBHs. Model data is analysed using a pipeline 
analogous to that used to reduce real cluster data. 

Section [2] gives a description of the models used in 
this work including the initial setup of the models and an 
overview of the evolution. A detailed look at the internal 
structure of the model clusters is then given in Section [3] 
along with a description of the attempt to analyse model 
data as real data. This is followed by a discussion in rela- 
tion to previous work and observations of Galactic GCs, and 
finally a summary of the main results. 



2 MODELS 

The focus of this work is a set of realistic A-body simula- 
tions that each starts with N — 100 000 objects - an object 
being either a star or a binary. Specifically, the starting mod- 
els contain: 100 000 single stars and no primordial binaries 
(labelled the K100-00 simulation); 95 000 single stars and 
5 000 binaries (K100-05); and, 90 000 single stars and 10 000 
binaries (K100-10). Masses for the stars are chosen from 
the initial mass function (IMF) of Kroupa, Tout & Gilmore 
(1993) between the limits of 0.1-50 M Q . Metallicity is set at 
Z — 0.001 for the stars. The initial positions and velocities 
are assigned according to a Plummer density profile (Plum- 
mer 1911; Aarseth, Henon & Wielen 1974) in virial equilib- 
rium. A scale length of 8.5 pc is set for each simulation - this 



is to comply with the tidal radius set by the external tidal 
field (see below). In actual fact the results from two sim- 
ulations starting with 100 000 single stars will be utilised. 
These simulations are identical in all respects except for the 
random number seed used to generate the starting masses, 
positions and velocities. These will be known as KlOO-OOa 
and KlOO-OOb. See Table [T]for a list of the simulations used 
in this work. 

The model clusters are evolved using the NB0DY4 code 
(Aarseth 1999, 2003). This includes algorithms for stellar 
and binary evolution as described in Hurley et al. (2001). 
Simulations are performed using 32-chip GRAPE-6 boards 
(Makino 2002) located at the American Museum of Natural 
History. Each simulation took approximately six months to 
complete on a dedicated GRAPE-6 board. 

To account for the tidal field of the Galaxy each cluster 
is placed on a circular orbit at a distance of 8.5 kpc from the 
Galactic centre with an orbital speed of 220 kms -1 . This 
is commonly referred to as a Standard Galactic tide (see 
Giersz & Heggie 1997 for a full description). For the model 
clusters in this work, which each have a starting mass of 
M ~ 50 000 Mq , this gives an initial tidal radius of about 
50 pc. With the length-scale given above the clusters are 
close to filling their tidal radii at birth, noting that the po- 
sition of the outermost star will vary from model to model 
as positions are drawn at random from a distribution. 

Each cluster was evolved to a minimum age of 16 Gyr. 
This ensured that the core-collapse phase of evolution was 
completed and that models of comparable age to GCs were 
available for analysis. In fact, for model KlOO-OOb it is not 
necessarily true that core-collapse was reached. For reasons 
that will become evident in Section [3] this model did not 
show a deep minimum in core-radius prior to its termination 
at 16 Gyr whereas the other three models did show such a 
minimum between 15 — 16 Gyr. For interest sake the K100-05 
simulation was allowed to proceed to 20 Gyr. After 16 Gyr of 
evolution the model clusters had been reduced to N ~ 22 000 
and, in terms of mass, approximately 80% of the cluster had 
been lost over that period. The tidal radius at 16 Gyr was 
about 30 pc. 

The evolution of the K100-05 model is shown in Figure[T] 
in terms of the number of half-mass relaxation times that 
have elapsed. This is done using both the initial half-mass 
relaxation timescale (t r h.o = 1 400 Myr) and the timescale 
after 15 Gyr (t T u t i5 = 580 Myr). The difference between the 
two is significant and shows that one must be very careful 
using the observationally determined of a cluster to infer 
the dynamical age (this point will be returned to later). The 
relaxation time is calculated according to the standard ex- 
pression developed by Spitzer (1987: see eq. 1 of Baumgardt, 
Makino & Hut 2005). In reality t r h is an evolving quantity, 
generally decreasing with age, and this must be accounted 
for when calculating the true dynamical age (see the solid- 
line in Figure QJ. The evolution of t T h for the K 100-00 and 
K100-10 simulations differs from the K100-05 simulation by 
no more than a few per cent across the evolution. 

3 RESULTS 

The evolution of r c /rh for the models starting with 0%, 5% 
and 10% binary frequency is shown in Figure [2] Here r c is 
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the density-weighted core-radius (Casertano & Hut 1985) 
commonly used in TV-body simulations and rh is the half- 
mass radius. These are not directly comparable to observed 
quantities. 

Initially r c /rh increases for all models. This is because 
the early phase corresponding to rapid mass-loss from mas- 
sive stars leads to an overall expansion and the effect is 
greater at smaller radii. Subsequent evolution has r c /rh gen- 
erally decreasing as it is dominated by the contracting core 
- rh continues to expand until about 4 Gyr. After that time 
the half-mass radius begins to feel the effect of the decreas- 
ing tidal radius and gradually decreases from that point on. 

The evolution of r c /rh is similar for all models at all 
times. At the 16 Gyr end-point there is some distinction 
between the models with and without primordial binaries: 
rc/rh ~ 0.07 in the former and ~ 0.02 in the latter. How- 
ever, the data in Figure [2] have been smoothed considerably 
using a moving 500 Myr window and 100 Myr increments. 
In reality the noise in the data would preclude drawing any 
inference regarding the primordial binary content of a clus- 
ter based on its r c /rh measurement. This is not to say that 
a systematic difference in core radius would not develop if 
the models were allowed to evolve well in to the post-core- 
collapse regime. 

In terms of comparing to real data the results of Fig- 
ure [2] are not particularly useful. What is needed is a pro- 
cedure that analyses the model data in the same way as is 
done for observations of clusters. In this way a meaningful 
rc/rh ratio can be extracted. The TV-body stellar evolution 
algorithm (Hurley, Pols & Tout 2000) provides the mass, 
luminosity and effective temperature of each model star. 
Using the model atmosphere data of Kurucz (1992), sup- 
plemented by Bergeron, Wesemael & Beauchamp (1995) for 
white dwarfs, these are then converted to broadband UVBRI 
colours. It is then relatively simple to calculate the half-light 
radius, r^i as the radius which encompasses the inner half of 
the total light of the cluster. This is a projected radius calcu- 
lated using a 2-dimensional projection of the 3-dimensional 
positions of the model stars. Finding the observational core- 
radius, which will be labelled r Ci i, requires analysis of the 
cluster surface brightness profile (SBP). For this it is possible 
to use the software described by Mackey & Gilmore (2003) 
in their work on the star clusters of the Large Magellanic 
Clouds. Each TV-body snapshot is taken in turn and used 
to construct a two-dimensional projected SBP. Stars more 
than two magnitudes brighter than the main-sequence turn- 
off and low-mass stars with Mv > 10 are excluded - this 
mimics the observational process of avoiding bright stars, 
which may saturate, and faint stars which may be incom- 
plete in number. Note that the projection is taken along the 
Y-axis and a choice is made to focus on the V magnitude. 
Neither of these choices affects the results to any signifi- 
cant degree. Next a three-parameter Elson, Fall & Freeman 
(EFF: 1987) model is fitted to the cluster SBP to deter- 
mine r c ,i (Mackey & Gilmore 2003). A similar approach was 
taken by Heggie et al. (2006) although the fit was made to 
the three-dimensional density profile and a fourth parameter 
was added in order to fit the central cusp for models with a 
central BH. 

As an example the SBP and EFF model fit for the K100- 
05 simulation at 15 Gyr is shown in Figure^. The resulting 
core radius is r c ,i = 0.99 pc. For comparison Figure|3]3 shows 



the projected surface density profile of the same stars along 
with the best fitting King model (King 1966). This gives 
r c> \ — 0.95 pc in good agreement. The corresponding TV- 
body core radius for the model cluster is 0.4 pc. Values of 
r c , rh, r c ,l (from EFF) and rh,i at 15 Gyr for each simulation 
are given in Table [T] 

Figure [4] demonstrates the relationship between the TV- 
body and observationally determined radii for the K100- 

05 simulation as it evolves. For the most part rh — 2rh.i 
in agreement with Baumgardt, Makino & Hut (2005). It 
is important to emphasize that r^,i is derived from a 2- 
dimensional projection of the TV-body data whereas rh is 
based on the original 3-dimensional data. Simply calculat- 
ing rh from a 2-dimensional projection gives a reduction of 
about 25%, as expected (see Fleck et al. 2006), and using 
the stellar light gives a further reduction. It can also be 
seen from Figure [3] that for the first ~ 7 Gyr r c is a good 
approximation to r c ,i. However, as the cluster becomes dy- 
namically old (t > 5t r h,o) this approximation is no longer 
valid. During core-collapse the central density increases and 
the value of r c computed from the density-weighted proce- 
dure decreases (as witnessed in Figure [2). At the same time 
the remnant fraction in the core is increasing (Baumgardt & 
Makino 2003) which flattens the profile of the visible stars 
and causes r c ,i to be greater than r c . 

After repeating the SBP-fitting process for the full set 
of simulations Figure [2] is repeated but now using r c ,i and 
rh,i. The result is shown in Figure [S] This ratio, (r c /rh)[, 
can be compared to observational data. It is clearly evident 
that the ratio is higher than previously reported - at 15 Gyr 
rc/rh ~ 0.3 regardless of binary content and without invok- 
ing an IMBH. Also plotted in Figure[5]is a fourth simulation, 
KlOO-OOb. The setup for this model was identical to that of 
KlOO-OOa except for the seed of the random number genera- 
tor. However, unlike KlOO-OOa this alternate model formed 
a BH-BH binary in the core after 4 Gyr of evolution. The 
BH masses are 24 and 25M0 and the binary formed in a 3- 
body interaction with an initial period of 19 000 d. At 16 Gyr 
it was still present in the core with a period of 195 d. The 
energy generated in 3-body encounters between this binary 
and stars in the core acts to 'puff-up' the core and inflate the 
core radius. This is analogous to what Baumgardt, Makino 

6 Hut (2005) find when an IMBH is present in the core. The 
KlOO-OOb simulation maintains r c /rh ~ 0.6 — 0.7 throughout 
the evolution and is clearly distinct from the other models 
from about 11 Gyr onwards. 

For comparison, the first long-lived binary in the KlOO- 
OOa simulation also formed at about 4 Gyr - comprised of a 
white dwarf and a helium star - but this simulation did not 
form a BH-BH binary at any point. The reason for this is re- 
lated to the velocity kicks given to supernovae remnants. For 
the models in this work, when a neutron star or BH is born a 
velocity kick chosen at random from a uniform distribution 
between — 100 kms is applied. This leads to retention 
fractions of 15-20% which is in line with the suggestions of 
Pfahl, Rappaport & Podsiadlowski (2002) for GCs. The ini- 
tial KlOO-OOa model contained 39 main-sequence stars with 
mass in excess of 2OM0, i.e. stars that would evolve to form 
BHs. However, only five BHs were retained in the model 
cluster after birth and only one of these BHs had a mass in 
excess of 2OM . The KlOO-OOb model started with 42 mas- 
sive main-sequence stars and had eight retained BHs, three 
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of which where more massive than 20Mq. So the velocity 
kick process, which itself is uncertain, and the related small 
number statistics of BH numbers, are certainly playing a 
role in determining the evolution histories of the models. 

In Figure [6]the projected surface density profiles of the 
KlOO-OOa, KlOO-OOb and K100-05 simulations at 15Gyr are 
compared. Profiles are constructed using radial bins of 500 
stars each and all stars are included. The profile of the 
K100-10 model is similar to that of the K100-05 model. 
Comparison of the KlOO-OOa and K100-05 profiles shows 
the expected result in that the single star model is more 
centrally condensed and would return a smaller core radius 
from King model fitting. By contrast the KlOO-OOb profile 
is much flatter. Thus the behaviour seen for iV-body mod- 
els with a central IMBH (Baumgardt, Makino & Hut 2005) 
can be replicated by the presence of a central stellar mass 
BH-BH binary. 

Finally, the KlOO-OOa and K100-10 simulations are used 
to look at the effect of primordial binaries on the distribution 
of remnants in evolved clusters. Baumgardt et al. (2003a) 
showed that the density profile of remnants (white dwarfs, 
neutron stars and stellar-mass BHs) rises more strongly in 
the centre of a cluster than the profile of luminous, or observ- 
able, stars. Thus the mass-to-light ratio rises naturally to- 
wards the centre of a cluster without the need for an IMBH. 
This point was shown by Baumgardt et al. (2003a) to be im- 
portant when interpreting the observed velocity dispersion 
profile of M15 which had been used to infer the presence 
of an IMBH in the core (Gerssen et al. 2002). The Baum- 
gardt et al. (2003a) models did not include primordial bina- 
ries. Thus the KlOO-OOa model in this work can be expected 
to show similar behaviour. Figure [T^, shows the projected 
density profiles of remnant stars and luminous stars (main- 
sequence stars with My < 10 and giants) at 15 Gyr. Indeed 
the remnant profile of model KlOO-OOa rises more steeply 
towards the centre. For this model at 15 Gyr remnants com- 
prised 40% of the cluster mass but only 1% of this was in the 
form of neutron stars and BHs. Note that to probe deeper 
into the centre of the model clusters 100 stars per bin has 
been used in Figure [7] which explains why the profiles are 
more erratic than those of Figures[3]D and|6] In Figure[7p this 
exercise is repeated for the K100-10 model - the presence of 
10% primordial binaries has erased any difference between 
the profiles. This is because binaries present a population 
of comparable average mass to the remnants and therefore 
segregate towards the centre on a similar timescale. 



4 DISCUSSION 

This work has gone some way to fulfilling a need identified 
by Trenti (2006) - taking realistic JV-body models and an- 
alyzing the snapshot data as if it were data acquired by a 
telescope. While these models are comparable in size to GCs 
at the lower end of the GC mass function (e.g. Gnedin & Os- 
triker 1997) they should not be taken as directly applicable 
to GCs. The results presented are mainly for comparison to 
other iV-body models - they provide an excellent companion 
to the models of Baumgardt & Makino (2003) and Baum- 
gardt, Makino & Ebisuzaki (2004) and a step forward in 
particle number compared to Trenti, Heggie & Hut (2007). 
Good agreement is also found with the Monte Carlo models 



of cluster evolution performed by Fregeau & Rasio (2007). 
Across a series of models starting with 100 000 objects (stars 
and binaries) these authors report r c /r^ values in the range 
0.05 — 0.1 with little to no dependence on the initial cluster 
profile or binary fraction. This is the same range shown for 
the iV-body models in Figure[5] noting that the Monte Carlo 
r c /rh is calculated using the traditional iV-body method. 
Significantly, Fregeau & Rasio (2007) do not see any no- 
ticeable change in r c /r-^ when they move to models starting 
with 300 000 objects. 

Trenti (2006) looked at globular cluster data from the 
catalogue of Harris (19983) to examine the distribution of 
r c /rh ratios. This involved carefully selecting a sample of 57 
GCs with the main determinant being that the measured 
half-mass relaxation timescale be less than 10 9 yr. This was 
to ensure that the GCs in the sample were all dynamically 
old - at least 10 half-mass relaxation times old based on a 
conservative age estimate of 10 Gyr. The motivation for do- 
ing this was based on the demonstration by Trenti (2006) 
that r c /rh can only be used to distinguish between clus- 
ters with differing initial content (single stars, primordial 
binaries and IMBHs) after at least 10 half-mass relaxation 
times have elapsed. However, this model result was based on 
the use of the initial half-mass relaxation timescale (t r h,o) 
whereas the observed clusters only give the current t T h- As 
has been shown in this work, this can be expected to be 
at least a factor of two less than i r h,o- Another consider- 
ation is that i r h calculated from models is using the three- 
dimensional half-mass radius while the value of f r h quoted in 
the Harris (1996) catalogue is based on the two-dimensional 
half-light radius. This can also lead to an overestimate of 
the true dynamical age of a clustei0. Thus caution is urged 
when comparing dynamical ages of model and real clusters. 

In fact, Trenti (2006) also considered a refined sample 
based on a criterion of i r h < 0.5 x 10 9 yr which is more 
appropriate in terms of ensuring the clusters are dynami- 
cally old. This led to a sample of 25 Galactic GCs. Of these 
there are 10 with r c /rh > 0.2 which is the condition used by 
Trenti (2006) to infer the possible presence of an IMBH. The 
larger sample considered by Fregeau et al. (2003: see their 
Fig. 17) shows that in general a Galactic GC is as likely to 
have r c /r\ > 0.2 than not. An important point here is that 
Fregeau et al. (2003) and Trenti (2006) are using r c /rh from 
models (corresponding to Figure [2} and this will underesti- 
mate the true ratio when compared to observations which 
use r Cj i/rh.i (as given in Figure |5). The results presented 
here show that observed ratios up to at least 0.4 do not re- 
quire an IMBH for explanation. The median r c /rh for the 
Trenti (2006) sample is 0.28 which is well matched by the N- 
body models. As a result it is suggested, conservatively, that 
rc/rh > 0.5 be used to distinguish clusters which require the 
presence of something out of the ordinary to explain their 
inner structure. In the Harris (1996) database there are six 
such clusters (with t T ^ < 0.5 x 10 9 yr). They are E3, Terzan 
3, NGC 6366, Pal 6, NGC 6535 and Pal8. 

1 The updated version of this catalogue is available on-line 
at either http://physwww.mcmaster.ca/%7Eharris/mwgc.dat or 
|http: / / coihue.rutgers.edu~andresj / gccat.html | 
z This point came to mind after noting the conversion applied 
in Baumgardt, Makino & Hut (2005) when calculating the relax- 
ation timescale for Galactic GCs 
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The iV-body models have shown that one explanation 
for a cluster observed to have a large r c /rh ratio is the pres- 
ence of a stellar mass BH-BH binary. The action of this BH- 
BH binary causes r c /rh to diverge from that found in models 
without such a binary. This divergent behaviour occurs after 
about 11 Gyr, in terms of model age (see Figure [SJ), which 
in dynamical terms equates to approximately six half-mass 
relaxation times (see Figure [T} . So it is possible to differ- 
entiate between models with and without a BH-BH binary 
before the completion of core-collapse evolution. 

The findings relating to the BH-BH binary model 
(KlOO-OOb) occurred very much by chance as this study in 
no way set out to create a model that would form such a 
binary. What it does demonstrate is that the random meet- 
ing of stars in a cluster (model or real) , and additionally the 
randomness introduced by velocity kicks given to supernova 
remnants, can have severe implications for the long-term 
structure and observed nature of a cluster. This particular 
model also showed a flattened density profile compared to 
the models that did not form a long-lived central BH-BH bi- 
nary. Baumgardt, Makino & Hut (2005) found that models 
with an IMBH also gave flattened profiles. Comparison to 
the observed surface brightness profiles of 37 Galactic GCs 
presented by Noyola & Gebhardt (2006) lead to the sug- 
gestion of five clusters of interest in terms of detecting an 
IMBH. These clusters may also be of interest for finding a 
BH-BH binary. 

The results in this study have actually made it more 
difficult to explain clusters with low r c /?"h - approximately 
half of the Galactic GCs have ratios less than 0.2 (Harris 
1996; Fregeau et al. 2003; Trenti 2006) and these cannot be 
reached by the models (based on inspection of Figure [5}. 
However, there are two factors to note here. The first is 
that the results shown in Figure [S] are smoothed. Looking 
at the raw, non-smoothed, data there is much fluctuation 
and values below 0.2 do occur in the models at late times 
(excluding the KlOO-OOb model) - the average error in the 
smoothed (r c /r^) 1 as shown in Figure [5] is approximately 
±10%. So small values are certainly possible depending on 
when a cluster is 'observed'. A word of caution is required 
on this point as real GCs are, in most cases, richer than 
the models presented here and statistical fluctuations will 
be smaller. The second point is uncertainty in the r c ,i fit- 
ting process. To demonstrate this one can look at the sur- 
face density profile of model K100-05 at 15 Gyr, as shown 
in Figure [3^. The King model fit shown gives r c = 0.95 pc 
however, if the fitting process is biased to fit the inner 1 pc 
of the profile then values as low as r c = 0.7 pc are plausible. 
It is also true that many of the GCs with r c /r^ < 0.2 in 
the Harris (1996) catalogue are also flagged as core-collapse 
clusters and accurate measurement of r c ,i using SBP-fitting 
software can be difficult for clusters passing through this 
phase. The interested reader can look at the SBP library in 
Trager, King & Djorgovski (1995) along with the associated 
description of the fitting process for core-collapse clusters. 

One could also expect that reducing the strength of the 
tidal field would lead to a reduction in r c /rh through an 
increase in tv However, there is no evidence in the Galac- 
tic GC sample for a link between r c /rh and distance from 
the Galactic centre. Note also that the models of Baum- 
gardt, Makino & Hut (2005) that are of comparable size to 
those presented here, but with an IMBH in the core, actu- 



ally give smaller r c /rh values. One difference between the 
two sets of models is that the Baumgardt, Makino & Hut 
(2005) models are isolated and indeed they do show a larger 
half-light radius. Even after correcting for this, the Baum- 
gardt, Makino & Hut (2005) models with an IMBH would 
give comparable (not larger) r c /?~h values to the models in 
this work without an IMBH. So it is not clear from this 
comparison that clusters with an IMBH should necessarily 
show a larger r c /rh ratio, as suggested by the models de- 
scribed in Trenti (2006) and theoretical arguments (Heggie 
et al. 2006). It is interesting to note an apparent discrepancy 
between the r c /rh values reported from the IMBH models 
of Trenti et al. (2007) and those of Baumgardt, Makino & 
Hut (2005). Modelling time-dependent tidal fields may also 
be important in determining the actual r c /rh ratio, and the 
choice of initial conditions, such as the scale radius, may 
also play a role. Clearly there is more work to be done in 
this field before we can resolve the issue of which GCs may 
harbour an IMBH. 

On one hand this investigation is suggesting that the 
presence of IMBHs in GC cores is not so likely - interme- 
diate r c /rh values can be explained by models without an 
IMBH provided the correct comparison is made and higher 
r c /rh values may instead show the presence of a stellar-mass 
BH-BH binary. However, the models also provide an oppor- 
tunity to look at how the distribution of stars in an old 
cluster is affected by the presence of a sizeable primordial 
binary population. The model with 10% primordial binaries 
shows that the mass distribution follows the light distribu- 
tion throughout the cluster - the steeper density profile of 
remnant stars compared to bright stars seen in the centre 
of single-star models is not replicated. Therefore, observa- 
tions that infer a steepening mass-to-light ratio in the core 
of a globular cluster should not be dismissed as a possible 
IMBH indicator (see also Gebhardt, Rich & Ho 2005). This 
is provided we assume that globular clusters are born with 
a modest binary fraction (Hut et al. 1992). Exactly how a 
direct model of a GC with primordial binaries and an IMBH 
will behave is beyond the scope of this work. 



5 SUMMARY 

By treating model data as if it were observational data 
higher r c /rh values than previously reported have been re- 
vealed. This provides a good match to the majority of the 
Galactic GCs without the need for an IMBH. It has also 
been shown that factors such as the presence of a BH-BH 
binary (comprised of stellar mass BHs) in a cluster core can 
flatten the measured luminosity profile and inflate the mea- 
sured core-radius. None of this precludes the existence of 
IMBHs in GC cores. However, it does demonstrate that the 
r c /rh ratio cannot be used with any certainty to infer the 
dynamical history or content of a cluster core. 
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Ratios of Star Cluster Radii 



integrated 

t-Arh.0 

--- t/W„ 

K100-05 / 




t/Myr 

Figure 1. A comparison of methods used for calculating the 
dynamical age of a star cluster in terms of the number of half- 
mass relaxation times elapsed as a function of time. Shown is the 
age of the cluster simply divided by cither the initial half-mass 
relaxation timescalc (t h = 1 400 Myr: dotted line) or the half- 
mass relaxation timescale at an age of 15 Gyr (th o = 580 Myr: 
dashed line). These are compared to the more detailed method of 
accumulating, or integrating, the number of half-mass relaxation 
times elapsed as the model evolves (solid line). Data is from the 
simulation starting with 95 000 single stars and 5 000 binaries. 
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Figure 2. Evolution of the ratio of the core-radius, r c , to half-mass radius, r^, for models starting with 0, 5 and 10% binaries (see 
Table [T] for a description). The radii are calculated using standard TV-body methods and three-dimensional data (see text for details). 
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r/pc r/pc 

Figure 3. Demonstration of the fitting process used to determine the observational core radius, r c i, using the K100-05 model at an 
age of 15Gyr as an example. Shown are: a) the V magnitude surface brightness profile with the best fit Elson, Fall & Freeman (1987) 
model; and, b) the surface density profile with the best fit King (1966) model. In both cases the data are projected along the Y-axis and 
stars fainter than My = 10 or more than two magnitudes brighter than the main-sequence turn-off are excluded. 
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5000 10 4 1.5x10 



t/Myr 

Figure 4. Comparison of radii calculated using the standard .W-body method to those calculated from fitting to the simulated luminosity 
profiles. Data from the K100-05 simulation starting with 95 000 single stars and 5 000 binaries. 
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Figure 5. Evolution of the ratio of the core-radius, r c 1, to half-light radius, i, for models starting with 0, 5 and 10% binaries (see 
Table [T] for a description). The radii are calculated from the simulated luminosity profiles using two-dimensional projected data (see text 
for details). An additional model that started with 0% binaries but formed a BH-BH binary at 4Gyr is included (KlOO-OOb). 
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Figure 6. Projected surface density profiles at 15 Gyr for the KlOO-OOa, K100-05 and KlOO-Ob models. All stars arc included and the 
profiles are computed using radial bins of 500 stars. 
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Figure 7. Comparison of the projected density profiles for remnant (circles) and bright (crosses) stars in: a) the KlOO-OOa model at 
15Gyr; and, b) the K100-10 model at 15Gyr. 
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Table 1. Parameters of the simulations performed in this work. Column 1 gives the label assigned 
to each model while Columns 2 and 3 show the starting number of single stars and binaries, 
respectively. Columns 4-7 give the following radii (in pc) for the models at an age of 15 Gyr: the 
.W-body density-weighted core radius; the half-mass radius; the EFF-fitted core radius; and, the 
half-light radius. Note that r c i and i arc from two-dimensional projected data while r c and 
arc based on three-dimensional data. 



Label 


N B 




r c 


rh 


r c ,l 


r h,l 


KlOO-OOa 


100 000 





0.34 


4.89 


0.85 


2.34 


KlOO-OOb 


100 000 





1.27 


5.59 


1.88 


3.72 


K100-05 


95 000 


5 000 


0.40 


5.25 


0.99 


2.75 


K100-10 


90 000 


10 000 


0.48 


5.31 


0.86 


2.71 



